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DEVELOPMENT  OF  SIDEBANDS  IN  TAPERED  AND  IN  UNTAPERED 

FREE-ELECTRON  LASERS 


I.  Introduction 

In  a  f ree-electron  laser  (FEL)  the  synchrotron  oscillation  of 
electrons  trapped  in  the  ponderomotive  potential  well  may  couple  energy 
into  sideband  frequencies.  The  ensuing  instability  leads  to  the  modulation 
of  the  output  signal  and,  in  consequence,  to  an  increase  in  its  spectral 
width. 

The  growth  of  sideband  frequencies  has  been  the  subject  of  discussion 
in  a  number  of  papers.^  ^  In  the  work  presented  herein  this  process  is 
examined  by  means  of  a  time-dependent  code  in  an  assumed  axisymmetric 
geometry.  Both  tapered  and  untapered  wigglers  are  examined.  Two  regimes 
of  sideband  development  are  examined  in  detail.  In  one,  the  carrier 
amplitude  is  small  and  all  frequencies  within  the  linear  gain  spectrum 
develop  independently.  In  the  other,  the  initial  amplitude  of  the  carrier 
is  large  and  hence  coupled  to  the  sideband  modes  via  the  synchrotron 
oscillation  of  the  electrons.  Simplified  analytical  models  of  sideband 
growth  in  the  two  regimes  are  presented  as  an  aid  to  understanding  the 
important  features  of  the  simulations. 

It  should  be  emphasized  at  the  outset  that  the  numerical  work  reported 
herein  is  intended  to  contrast  the  development  of  sidebands  in  untapered 
and  in  tapeied  FEL  systems.  In  this  connection,  the  important  issue  of 
sidcnard  start-up  from  the  noise  spectrum  appropriate  to  the  electron  beam 
and  wiggler  parameters  employed  is  not  addressed. 


II.  Numerical  Model 


This  section  presents  the  equations  that  form  the  basis  for  the 

investigation  of  the  PEL  interaction. 

The  electrons  move  under  the  influence  of  the  ponderomotive  force  due 

to  the  beating  of  the  wiggler  field  with  the  optical  field.  The 

computations  are  simplified  by  employing  the  Gaussian-Laguerre  Source- 

8  9 

Dependent  Expansion  (SDE)  technique,  ’  thereby  minimizing  the  number  of 
Laguerre  polynomials  required  for  an  accurate  description  of  the  optical 
field . 

The  optical  field  is  taken  to  be  of  the  form 


(r,z,  t) 


=  2  a(r,z, t) 


exp 


c .  c . 


9 


2 

where  A^.  =  me  a^/ |e  j  is  the  radiation  vector  potential,  m  is  the  rest-mass 
of  an  electron,  |e|  is  the  magnitude  of  the  electronic  change,  c  is  the 
speed  of  light  in  vacuo ,  w  is  the  radian  frequency,  and  e^  is  the  unit 
vector  along  the  x  axis.  The  wiggler  field  is  assumed  to  be  plane- 
polarized,  of  amplitude  and  period  2n/k^: 


B  (z)  =  B  exp  (ik  z)e  +  c.c. 
-  w  2  w  w  -y 


where  transverse  variations  of  the  wiggler  field  are  neglected,  and  e^  is 

the  unit  vector  along  the  y  axis.  The  equations  of  motion  of  the  j-th 

2 

electron,  of  energy  Yj^c  ,  are  then  given  by 


I"!  Y  .  i  0‘ja  f  _ 
_  ^  _ V  _B 

fl  f  4  Y  . 


E 


a  L 
n  n 


(> 


(1-ia) 


c.c. 


(1) 

(2) 


2 


where 


(3) 


defines  the  resonant  (i.e.,  synchronous)  relativistic  factor,  a^ 

2 

=  je|B^/mc  is  the  normalized  vector  potential  of  the  wiggler,  r^  is  the 

radial  distance  of  the  i-th  electron,  =  (to/c  +  k  )  z.  -  oot  is  the 

J  ^  3 

relative  phase  with  z^  the  axial  location  of  the  j-th  electron,  and  fg  = 

J^(  E.)  is  the  usual  difference  of  Bessel  functions,  with  E  = 

2  2 

(a  /2)  /(1+a  /2).  The  results  to  be  presented  in  this  paper  pertain  to 

V  w 

the  case  where  only  a  is  tapered  and  the  period  of  the  wiggler  is  taken  to 


8 


be  constant. 

Following  Sprangle  et  al.,"^  the  envelope  of  the  radiation  field  is 
expanded  as  follows 


a(r,z,t)  =  Y,  a^(z,t)  2r^/r^^j  exp  [^-(l-ia)r^/r^^j  . 


Here,  a(z,t)  is  related  to  the  curvature  of  the  optical  wavefronts,  r^(z,t) 

2  2 

is  the  spot  size,  and  (2r  /r^  )  is  the  Laguerre  polynomial  of  degree  n. 
The  method  of  SDE  then  permits  a  complete  specification  of  the  optical 

g 

field  by  solving 


9  1 

0  2 

I  2c  a 

> 

lar 

^  azj 

r  =  -  -  r  c 

'  s  wr  s 

s 

‘ 

M 

2  2 

I  2c  (1+a  ) 

a  =  - ^ ^ 

2c  (n  - 

l8r 

3z  J 

'  / 

h- * 

1 

SC 

\ 

1 

(4) 


(5) 


''d  d  ^  * 

♦  r  ^  cA  a  ^  inBca  ^  +  i  (n+1)  B  ca  ,  -  icF  ,  (6) 

9z  n)  n  n-1  n+1  n 
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sideband  development  is  as  described  herein.  Finally,  in  all  the 
computations  the  initial  electron  distribution  is  taken  to  be 
monoenergetic,  the  radial  profile  of  the  electron  beam  is  taken  to  be 
parabolic,  and  betatron  motion  is  neglected. 


III.  Numerical  Results 


For  definiteness  the  parameters  for  the  computations  presented  herein 

correspond  to  those  of  the  Paladin  experiment  at  the  Lawrence  Livermore 

National  Laboratory , and  are  'isted  in  Table  I. 

In  the  linear  regime  the  maximum  growth  rate  obtains  at  the  resonance 

wavelength  =  10.34  pm.  In  what  follows,  where  appropriate,  a 

wavelength  X  will  be  denoted  in  terms  of  the  relative  shift  S  s  (X/X  )-l 
^  res 

from  the  resonant  value.  Note  that  the  wavelength  of  the  carrier  is  given 

by  2nc/(A).  Spectra  are  obtained  by  performing  a  spatial  Fourier  series 

analysis  of  the  optical  field  a(r=o,z,t)  along  the  mesh.  The  result, 

displayed  in  the  figures  as  Fourier  amplitude,  is  dimensionless. 

a )  Mult i -frequency,  small  input  power 

Simulations  of  a  pulse  of  radiation  extending  over  many 

ponderomot i ve  buckets  permit  the  development  of  frequencies  besides  that  of 

the  input  signal,  and  as  electrons  slip  relative  to  the  radiation  pulse  the 

sideband  instability  may  develop. 

Figure  1  shows  the  development  of  3  spectral  components  of  the 

_2 

optical  field  when  the  input  signal  is  at  10.6  pm  (6  =  2.51x10  )  and 

IMi)  w.  Of  course  there  are  many  other  spectral  components  besides  those 

shown  in  Fig.  1;  liowever,  the  curves  shown  do  indicate  the  general  trend  in 

I  lie  development  of  the  sideband  frequencies.  Note  in  particular  that  in 

_2 

the  e>:[)onen  t  ia  1  regime  the  growth  rate  of  the  component  at  S  =  0.77x10 

exceeds  that  of  tiie  main  signal.  This  result  is  discussed  in  Section  IV. 

h)  Mill  I  i  frcsniency^  untapered  magnet 

FigiiLe  7  .shows  the  development  of  the  carrier  for  the  case  where 

_2 

i'ld  MU  Ilf  10. o  (jtn  radial  if)n  (6  =  2.51x10  )  is  injected  into  an  untapered 

nap, net.  Again,  there  ate  many  inodes  in  the  spectrum  that  grow  along  the 
letigtii  of  the  /iKgler.  I'he  dashed  curve  in  Fig.  2  indicates  the  maximum 
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amplitude,  or  the  envelope,  of  the  rest  of  spectrum  as  a  function  of  z. 

Fig.  2  also  shows  in  detail  the  evolution  of  one  of  the  fastest  growing 

modes,  at  5  =  5.71x10  ,  indicating  the  trend  in  the  development  of  the 

instability.  At  the  end  of  the  wiggler  (25  m)  the  amplitude  of  the 

sidebands  is  large  enough  to  spatially  modulate  the  optical  field  by  about 

30°J.  Figure  3,  which  shews  the  phase  of  the  carrier  along  the  wiggler, 

will  be  discussed  in  connection  with  the  analytical  model  for  the 

instability  in  Section  IV.  At  the  wiggler  exit  the  electron  beam 

distribution  function  has  a  clear  multi-stream  character  as  indicated  in 

Fig.  4.  The  spectrum  of  the  optical  field  at  this  point  consists 

principally  of  the  carrier  with  a  group  of  Stokes  and  anti-Stokes  modes  on 

oitiiei  side,  Fig.  5.  The  approximately  symmetrical  form  of  the  spectrum, 

■-'iurh  is  a  reflection  of  comparable  growth  rates  for  modes  symmetrically 

disposed  with  respect  to  the  carrier,  is  discussed  in  Section  IV. 

r  )  Mu  1 1 i - 1 requency ,  tapered  magnet 

As  is  well  known,  for  practical  purposes  the  magnet  employed  in  an 

M’[  device  must  be  tapered  so  as  to  enhance  its  efficiency  and  extraction. 

f)  shows  tiie  development  of  the  carrier  and  of  the  maximum  amplitude 

v  -ill  the  sidebands  through  a  device  where  the  normalized  vector  potential 

u  ,  |c|b.  ■^me'^'k^  is  tapered  as  shown  in  Fig.  7.  The  form  of  the  tapering 

.f?d  in  the  computations  is  obtained  by  simply  prescribing  a  constant 

I  of  decrease  of  energy  for  a  synchronous  electron,  at  an  assigned 

1  id  i  u  ; .  Fi  om  F,(] .  (3)  vitli  dy  /dz  =  constant,  one  obtains  a  (z).  Comparing 

r  w 

'  ,  anri  fi  if  is  apparent  that  in  the  tapered  device  the  sidebands 

■  i!  liurt.  at  a  le’/oJ  substantially  below  tliat  of  the  carrier.  As  expected, 

•  i.i  wf)ii<al  field  is  observed  to  be  only  slightly  modulated  in  space.  The 
‘  !•  t-'iii  r|  1  : 1 1  i  bu  t  i  on  ftni'-tion  (not  shown  hete)  consists  principally  of  two 


IV.  Analysis  of  Results 


,*■ 


An  understanding  of  Fig.  1  may  be  obtained  by  performing  a  single- 
bucket  linear  stability  analysis  of  Eqs.  (l)-(7).  The  presentation  is 


limited  to  the  fundamental  optical  mode,  and  a  monochromatic  electron  beam 


of  energy  per  electron. 


Defining 


tj/.  =  ^  t , 

J  J  o 


r.  -  y./y  , 
J  jo’ 


=  ia  exp  (iii/^t), 


the  equilibrium  corresponds  to  A  =  0,  =  1,  Zexp(-inij/j )  =  0  (n  =  1,  2,...) 


Perturbing  Eqs.  (l)-(7),  defining  collective  variables  as  in  Ref.  12, 
and  assuming  a  temporal  dependence  of  the  form  exp  (-iTt),  the  following 
disi)ersion  relation  is  obtained 


-  2c^(l-ia)/wr^^ 


N 


T  -  2  ck 


oj 


(®w^b) 


Nr 


(9  ^ 


0  , 


K I e  a  2  E  (1  -  exp  (-2  r^^/r^^).  Perturbing  B  [defined  by 


J 


1>I .  (H)|,  and  numerically  solving  Eqs.  (4)  and  (5)  along  with  the  cubic  in 
T,  niK  n!)i,ains  the  giowth  rate,  efficiency,  spot  size,  and  a  in  the 
<■  xp'iiK  11 1  i  a  1  regime  for  any  given  angular  frequency  w. 

(iov,  as  is  well  known  the  FEL  interaction  lias  the  rather  important 
I-:  r^iK'i  t  y  Hia  f  tlie  optical  field  tends  to  be  guided  by  the  electron  beam. 

Ill  an  am[;lifier  operating  in  the  exponential  regime,  it  is  found  that. 


^  *1 


irrespective  of  the  initial  spot  size,  the  radiation  beam  asymptotes  to  a 

unique  spot  size  r^  and  wavefront  curvature  (~  a”^)  (Rets.  8  and  9). 

Figure  8  shows  the  growth  rate  T,  efficiency  ti,  matched  spot  size  r^  and 

matched  a  in  the  case  of  a  small  input  power  as  a  function  of  5  = 

(X/\^^^)-l,  where  X  =  2nc/(»).  The  crosses  are  the  results  of  single-bucket 

simulations.  The  curves  are  obtained  from  the  linear  stability  analysis  of 

the  preceeding  paragraph.  It  is  seen  that  the  agreement  is  quite  good. 

An  important  feature  of  Fig.  1  may  now  be  understood  with  reference  to 

Fig.  8.  In  the  small  signal  regime  -  and  therefore  prior  to  particle 

trapping  -  Eqs.  (.1)  -  (7)  may  be  linearized  to  show  that  there  is  no 

coupling  between  the  various  spectral  components.  In  other  words,  the 

development  of  the  spectral  components  proceeds  independently  and  at  a  rate 

approximately  equal  to  that  indicated  in  Fig.  8.  Referring  to  Fig.  1,  it 

_2 

is  thus  seen  that  the  larger  growth  rate  of  the  sideband  at  5  =  0.77.xl0 

_2 

as  compared  to  the  carrier  at  5  =  2.51.xl0  is  consistent  with  Fig.  8. 

It  is  also  possible  to  set  up  and  analyze  a  simple  model  of  the  FEL 
interaction  so  as  to  obtain  an  understanding  of  the  gross  features  of  the 
cideiiand  in^taliility  ensuing  from  the  synchrotron  oscillation  of  the 
ele<. trolls  in  a  large-amplitude  carrier  wave. 

Neglecting  diffraction  and  considering  the  fundamental  mode  of  the 
opiical  field  only,  writing  Yj  =  Yj.  +  Ayj  ( t ) ,  iJ/j  =  +  Aiiij(t),  a  =  | 

•  ;::p  the  equilibrium  is  described  by 


'I/.  ^  j>  .  ar.  /r  ^  2n .  JI 

)  o  ]  s  j 


(nj 


is  an  integer) 


A'i/.  -  -  9.“^  Aii/.  ,  (9) 

fit'  ’ 


■  finr,  t  an  t  , 


10 


c  3  ) 

♦  <“>  .  X.  ( 

2c)  y.  B  w 

Ut 

3zj 

.r  j  .  Y  j  ’ 

s  j  r 

where  Che  synchrotron  frequency  is  given  by 


(10) 


Q. 

J 


ck.  wa  fn 

w  V  B  I  (o) 
n  I  S 


,1/2 


e. 

J 


?  2 

£.  =  exp  (-r.  /r  ),  and  y  is  defined  by  Eq.  (3).  Note  that  Eq.  (10) 

J  J  2  ^ 

describes  a  linear  increase  in  the  equilibrium  phase.  Inserting  numerical 
values  corresponding  to  the  wiggler  entrance,  Eq.  (10)  yields  a  rate  of 
inciease  that  is  about  a  factor  of  4  greater  than  that  observed  in  Fig.  3 
(beyond  9  m).  Due  to  diffraction,  however,  the  spot  size  r^  increases 
along  the  wiggler.  Thus,  inserting  an  average  value  for  r^  (between  9  m 
and  25  m)  Che  rate  obtained  from  Eq.  (10)  is  within  several  tens  of  percent 
of  the  numerical  value. 

Next,  perturbing  the  equilibrium  state  and  assuming  a  dependence  of 
the  form  f^^^  --  exp  [-i(2.(A)t  -  Ak  z)]  for  the  perturbation,  one  obtains  the 
tol lowing  dispersion  relation  (Appendix  A) 


1 

X 


[leck.^  '  ^'b'] 

(.v-o' 


_6 _ 

X(X-C)  ' 


(11) 


f  ^  7 

L'.)  1  ^  a  ^/2  /2Q„y 

V  ;  B  r 


,  C  =  cAk 


(l  .  a,//2)  /20j- 


a  t 


,1  <'))  , 

Y, 


1  ^  ^  H': 

I,  w  s 


I  J  /N|  2 

)  i  j 


1 1 


=  E  Q  /N 


T^ig  dispersion  relation  is  obtained  in  the  diagonal  approximation,  wherein 


terms  of  the  form  f^^^  ,  f^^^  Ayj  are  neglected.  The  dispersion 


relation  in  Eq.  (11)  is  similar  to  that  given  in  Ref.  6.  The  last  term  on 


the  right-hand  side  of  Eq.  (11),  which  is  due  to  terms  proportional  to 


y.'  and  was  neglected  in  Ref.  6,  increases  the  growth  rate  by  about  10^. 


In  reference  to  the  electron  distribution  function  shown  in  Fig.  4,  it  is 


interesting  to  note  the  similarity  of  Eq.  (11)  to  the  dispersion  relation 


for  a  multi-stream  system.  Neglecting  the  third  term  on  the  right-hand 


side  of  Eq.  (11),  it  is  seen  that  -  as  discussed  in  Ref.  6  -  the  dispersion 


relation  has  a  symmetrical  form  about  the  carrier  (Ak.  =  0),  with  identical 


growth  rate  for  (-Aw,  -Ak)  and  (Aw,  Ak) .  The  spectrum  in  Fig.  5  is  in 


approximate  agreement  with  this  general  feature  of  Eq.  (11). 


Inserting  numerical  values  into  Eq .  (11),  one  finds  that  the  range  of 


unstable  modes  encompasses  that  obtained  in  the  simulations.  The  maximum 


growth  rate  obtained  is  about  a  factor  of  5  larger  than  the  average  value 


observed  for  one  of  the  fastest  growing  modes  in  the  simulations  [cf.  Fig. 


That  the  growth  rate  deduced  from  Eq.  (11)  should  exceed  that  observed 


in  the  simulations  and  the  discrepancy  in  the  rate  of  increase  of  the  phase 


of  the  carrier  are  to  be  expected  for  two  reasons.  First,  diffraction  of 


tlie  optical  field  [neglected  in  deriving  Eq.  (11)  |  is  bound  to  reduce 


i-/th  latec.  Second,  in  the  simulations  the  electrons  are  distributed 


thioughout  the  ponderomo t i ve  bucket  with  synchrotron  frequencies  ranging 


tiom  Q.  down  to  zero,  whereas  the  analysis  leading  to  Eq .  (11)  assumes  all 


ctioiis  t  f)  i)ounce  at  the  b(jttom  of  the  ponderomo  t  ive  wells,  at  the 


laige^'t  syuf  h  L 1)  I  1  frequency  (cf.  Eq  .  (9)]. 


Turning  next  to  the  case  of  the  tapered  magnet,  the  question  arises  as 

to  why  the  growth  rate  of  the  sidebands  is  about  a  factor  of  6  smaller  than 

that  in  the  untapered  magnet.  [cf.  Figs.  2  and  6].  As  is  well-known,  upon 

tapering,  the  electrons  separate  into  roughly  two  groups.  For  the  tapering 

employed,  the  decelerating  group  -  which  is  responsible  for  stimulated 

emission  -  comprises  25-30%  of  the  total  number  of  electrons.  However,  the 

sideband  growth  rate  for  such  a  fraction  of  the  number  of  electrons  is 

still  large  compared  to  the  observed  value.  There  appears  to  be  additional 

reasons  for  the  small  growth  rate  of  the  sidebands.  In  the  context  of  a 

high-extraction  FEL  with  a  tapered  magnet  it  is  generally  assumed  that  the 

ciiange  of  parameters  (a_  and/or  k  )  is  slow  enough  so  that  the  action 

V  y 

I.  H  rdii(.  (Y'-y  )  is  an  adiabatic  invariant.^  For  the  case  considered 

J  J  J 

here,  however,  it  turns  out  that  dY^/dt  =  ^^(Yj-Yj.),  implying  that  the 
action  is  not  an  invariant  (Appendix  B).  Additionally,  the  Yj  of 
trapped  electrons  decrease  more  or  less  monotonically  with  time,  modifying 
rhe  "equilibrium"  electron  distribution  on  the  same  time  scale  as  that  of 
tlie  synchrotron  motion.  Physically,  as  the  radiation  pulse  slips  relative 
fo  rhe  electrons,  it  is  modulated  not  only  by  the  synchrotron  motion  but 
also  1)7  rhe  temporal  variation  of  y^,-  This  has  a  detuning  effect  and 
renders  the  sidebands  quasi -stable.  Indeed,  by  increasing  the  rate  of 
tap^^’fing  and  therefore  the  rate  of  variation  of  y^.,  the  modulation  of  the 
r.Driral  pulse  is  reduced  further,  although  the  efficiency  is  also 
d  i  mi  n  i  .'.died  due  to  increased  particle  detrapping. 

Finailv.  it  is  important  to  remark  that  the  tapering  employed  here 
■oirs  "o  file  same  efficiency  for  a  ingle-ltucket .  single- f  requeue;,’  case  as 
f';t  .1  nairz-lujcke  r  ,  mul  t  i- frequency  case.  In  other  words,  for  the  gi’/en 

f')!m  for  a,,(c),  the  optimal  rate  of  tapering  naturall;,’  leads  to 
•;U'  '  i'S'-f've'i  rpriiicrion  in  the  growth  rate  of  sideband  modes. 


'  •  * 


V. 


Conclusion 


The  work  presented  here  contrasts  the  development  of  sideband 
frequencies  in  an  untapered  and  in  a  tapered  magnet.  For  the  tapering 
employed,  optimal  operation  (i.e.,  maximum  efficiency)  is  achieved  when  the 
electrons  trapped  in  the  ponderomotive  buckets  are  non-adiabatically 
decelerated,  leading  to  substantial  suppression  of  the  sideband  frequencies 
relative  to  the  untapered  system. 
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Appendix  A:  Multi-frequency  linear  stability  analysis 
In  this  appendix  some  of  the  details  leading  to  the  dispersion 


fr---- 


I 


relation  in  Eq.  (11)  are  presented. 

Upon  writing  y.  =  '^io 

k(o)  |a|  =  I  +  ja^^^l,  and  neglecting  terms  comprising 


=  r 


products  of  the  form  Ay.  A<|/.  etc.,  one  obtains,  at  first 

J  j  J  J 


order . 
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Mote  that  the  term  proportional  to  Yj  on  the  right-hand  side  of  Eq.  (Al) 


was  neglected  in  Ref.  6  where  a  similar  analysis  of  sideband  growth  was 
pi esen  t  ed . 

Defining  collective  variables 


Xj,  -  I  y  =  Es.’^y/^^/N,  E,,  =  Ee,‘^/N 
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g 

dt 


1_ 

dt 


yi 


w 


oi  a 

27" 

r 


+  E. 


» 


,  2ck 

d  ^  _ w 

dt  ^2  =  Y  ^2  ’ 
r 


etc. 

This  hierarchy  may  be  truncated  by  assuming  X2  =  Ze^X^/N,  E2  =  lEjE^/N, 
whence 


dt 


IlV 

2ck. 


Xi  +  Ej 


f 


where  Qg  is  defined  by  Eq.  (12). 

Assuming  perturbations  of  the  form  exp  [-i(Acot  -  Akz)],  one  obtains 
the  dispersion  relation  in  Eq.  (11),  where  the  third  term  on  the  right-hand 
side  is  due  to  the  term  proportional  to  on  the  right-hand  side  of 

Eq.  (Al). 
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Appendix  B:  Temporal  variation  of  adiabatic  invariant 
In  this  appendix,  the  law  of  variation  of  the  action  variable  for  the 
tapered  magnet  is  obtained. 

For  a  synchronous  particle,  the  energy  and  phase  evolve  according  to 


dT  ^r  = 


w  a  fo 
w  B 

2y. 


|a|£^sin(^g/^  +  (j)  +  > 


dt  ""r  =  0  . 

Defining  Yj  =  Yj.  +  Syj  ,  +  S'f'j ,  the  Hamiltonian  for  small  amplitude 

oscillations  is 


ck  «  w  a  f„  a 
w  ,  2  w  B  ' 

H  =  -  6y-  +  - '/ - 

Y  J  ^Y 

r  r 


/  2  2'\ 
e  cos  N/  +  <j)  +  ar  /r  . 

r  ^r  r  s>)  j 


Defining  the  synchrotron  frequency  of  such  a  particle, 

1/2 

Q  =  ck  0)  a  fo  la  I  c  cos  np  +  <|>  +  ar  /r  )  /y  , 
syn  L  w  w  B'  '  r  t  r  r  s/J  r  ’ 

following  Ref.  13,  the  action  (Ij)  angle  (w^)  variables  are  found  to  evolve 
according  to 


- 

dt  j 


I .  dy 

cos2w.  , 

Y  dt  1  ’ 

r 


d  1  '^^r 

jr  w.  =  +  = —  t;—  sin2w. 

dt  J  syn  2y^  dt  j 


These  formulae  are  valid  for  the  length  of  the  magnet  lying  between  6  m  and 

21  m  where  52  is  found  to  be  remarkably  constant  in  the  simulations, 
syn 

e  variation  of  y^^  is  the  only  factor  contributing  to  the  breaking  of 


Thus  th 


the  adiabatic  invariant  .  (Generalization  of  the  formulae  to  include  the 
variation  of  52^y^  iti  the  other  sections  of  the  magnet  is  straightforward.) 
The  change  in  the  action  variable  over  a  synchrotron  period  is  then  found 
to  be  given  by 
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It  is  important  to  notice  that  although  this  variation  is  relatively  small, 
it  is  nevertheless  much  faster  than  the  usual  case  for  an  adiabatic 
invariant  where  M  ->  0  exponentially  as  dYj./dt  -*  0. 
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